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Abstract 

The modified gravity with In R or R~'^ (In i?)™ terms which grow at 
smah curvature is discussed. It is shown that such a model which has 
well-defined newtonian limit may eliminate the need for dark energy 
and may provide the current cosmic acceleration. It is demonstrated 
that R^ terms are important not only for early time inflation but also 
to avoid the instabilities and the linear growth of the gravitational 
force. It is very interesting that the condition of no linear growth for 
gravitational force coincides with the one for scalar mass in the equiv- 
alent scalar-tensor theory to be very large. Thus, modified gravity 
with R^ term seems to be viable classical theory. 
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1. Introduction. The astrophysical data from high redshift surveys of 
type la supernovae [1] and from the anisotropy power spectrum of CMB 
[2] change our image of current universe which seems to be accelerating. 
The theoretical foundation which is used to construct such universe is also 
quickly evolving. In particular, the popular explanation of current universe 
acceleration is based on the dominance of some mysterious, exotic matter 
called dark energy. There is still no the satisfactory theoretical explanation 
for the origin of this exotic matter which should appear precisely at current 
epoch. 

Having in mind that more astrophysical data will be available soon, it 
seems the right time to search for alternative explanation of current cosmic 
speed-up. In recent papers [3, 4] (see also [5]) it has been suggested the 
gravitational alternative for dark energy. The key idea is to modify the 
Einstein action by l/i? term which dominates at low curvature. Remarkably, 
that such terms may be predicted by some compactifications of string/M- 
theory [6]. Unfortunately, it has been found [7, 8] that 1/R model contains 
unacceptable instabilities which does not appear in the Palatini version of the 
theory [9, 10, 11]. (It is known that in general case Palatini and standard 
metric versions are not the same, basically they lead to different physical 
systems [12]). Moreover, the scalar-tensor theory which seems to be non- 
realistic [13] in the standard formulation due to solar system observations 
[14] becomes viable in Palatini form. 

Nevertheless, more comphcated modification of Einstein gravity of the 
form B? -\- R -\- 1/R [15] predicts the unification of the early time infla- 
tion and late time cosmic acceleration in the standard , metric formulation. 
Moreover, the instability found in ref. [7] is significally suppressed by higher 
derivative(HD) R^ term or other HD terms like R^. The solar system test for 
equivalent scalar-tensor theory [13] may be passed because scalar has large 
mass induced again by HD terms. The consideration of such theory in Pala- 
tini form has been recently done [16]. It is shown that account of R^ term 
makes the theory viable also in Palatini form. 

In the present paper we continue the search for the realistic modified 
gravity which may provide the gravitational alternative for dark energy. As 
such a model, we suggest to account for the In it! terms in modified gravity. 
Such terms are basically induced by quantum effects in curved spacctime. 
Various versions of such modified gravity may eliminate the need for dark 
energy and may serve for the unification of the early time infiation and cosmic 
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acceleration. HD terms again suppress the instability and improve the solar 
system bounds so that the theory may be viable. The correction to the 
gravitational coupling constant may be small enough too. 
2. General formulation and simplified model. One may start from 
rather arbitrary function F{a,b^,c, g^i,) which depends on two scalars a, c, 
one vector and metric g^,^. The general starting action is: 

S^^J d^xy^F{R, d^R, OR, g^,) . (1) 

Here R is the scalar curvature. Introducing the auxiliary fields A and B, one 
may rewrite the action (1) as following: 

S^^J d'xV^ {B {R-A) + F{A, d,A, DA, g^,)} . (2) 

By the variation over B, A — R follows. Substituting it into (2), the action 
(1) can be reproduced. Making the variation with respect to A first, we 
obtain 

B = f{A) = {daF{a, b^, c, g^^) - {db^F{a, b^, c, g^^)) 

which may be solved with respect to A a.s A = G{B). In general, A is solved 
non- locally as a function of B. Eliminating A in (2) by using G{B), we obtain 

S^^j d'x^g{B{R-G{B)) + F{G{B),d^G{B),nG{B),g^,)} . (4) 

Instead of A, one may eliminate B and arrive at the equivalent action. 

The scalar field a may be defined a = — Ini? = In /(A). One can scale 
the metric by g^i, — > e'^g,j,v Then the action (4) can be rewritten as 

S = ^ j d''x^^^R-'^-g^'''d^ad,a-V{a)^ 

V{a) = G'(e-'^)e^-e2'^F(G'(e-'^),9^(G'(e-'^)), 

□ (G (e-'^)) + g^^d.ad, [G (e"'^)) , e^g,,) . (5) 

It is given in the Einstein frame. On the other hand, the (physical) action 
(4) is given in the Jordan frame. 
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As the simple example, we consider the following case: 



F{A)^A + aln^. (6) 
A* 



One finds 



e- = S = l + an( — ) , (7) 



which can be solved non- locally with respect to A. Then the complicated 
expression for the potential follows 



V{a) = e'^D-i ( -—-^ . ) - e^- 



°"' I □-(e-»-l) .' 



+aln I— — — ^ + g^^'d.ad, ( Q-^ ( , ^ 

We may consider more general example 

F{A) = A + a]n — + a'ln — + pA"" . (9) 
/J, /J, 

This model is very complicated. If we consider the case R is almost constant, 
the second term turns to the (cosmological) constant. The natural starting 
model looks as follows: 

F(A)=A + a'\nA, + PA"'. (10) 

Furthermore, m = 2 choice simplifies the model. The generalizations of 
this model will be considered at final section. The correspondent term 
at large curvature leads to well-known trace anomaly driven (Starobinsky) 
inflation. Assuming the scalar curvature is constant and the Ricci tensor 
is also covariantly constant, the equations of motion corresponding to the 
action (1) are:^ 

= 2F(R) - RF'iR) = f(R) = R + 2a' In ^ - a' . (11) 



^If the action contains only term, F{B) = /3i?^, Eq.(ll) is trivial, which means 
the arbitrary constant curvature space is a solution. Adding the Einstein-Hilbert term, 
F{R) = R+pB?, the only solution of (11) is i? = 0. Then if we start with large R solution, 
which may correspond to the inflation, due to the Einstein-Hilbert term, R decreases slowly 
and goes to zero. The inflation will be stopped. 
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If OL > 0, /(-R) is monotonically increasing function and we have 
liniij^o /(-R) = —oo and liniij^+oo /(-R) = +00. There is one and only 
one solution of (11). This solution may correspond to the inflation. On 

the other hand, if a' < 0. limR-^o f{R) = 00 and lim/j^_(_oo = +00. Since 
f'{R) = 1 + ^, the minimum of /(-R), where f'{R) = 0, is given by 
R = -2q'. If /(-2a') > 0, there is no solution of (11). If /(-2a') = 0, 
there is only one solution and if f{—2a') < 0, there are two solutions''. Since 

/(—2a') = —2a' ^1 — In^^j, there are two solutions if — ^ > e. Since 
the square root of the curvature R corresponds to the rate of the expansion 
of the universe, the larger solution in two solutions might correspond to the 
inflation in the early universe and the smaller one to the present accelerating 
universe. 

For m = 2 case, we have 

e-- = l + ^ + 2/3A, (12) 
which can be solved with respect to A: 



- (1 - e-'^) ±\ (1- e-"? - 8/?a' 

■ 

In the branch that the r.h.s. in (12) is negative, one may choose 

e-'^^-(^l + ^ + 2pA^ ' 

Then instead of (5), we obtain 

S = ^Jd^x^i^-R+^g>'''d^ad,a-V{a)'^ 

V{a) = -G (e-'^) e^ - e^'^F (G (e-^)) . (15) 

As the sign in front of the scalar curvature R is changed, this seems to be 
the indication to the anti-gravity. Of course, the anti-gravity should not be 



^It is interesting that only In R may lead to early time inflation and current acceleration. 
The AdS cosmologies [17] may be considered as well. 
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real but apparent since the physical theory in the real spacetime should be 
given by (1) with (10) (m = 2 case). 
The potential is 

a' (l - In 4) + 
ViA) = V(A) = — ^ ^ . (16) 

Then in terms of cr, V{a) can be expressed as 
V{a) = e^'^ J a' I 1 - In 



- (1 - e-'^) ± - e"^)^ - W 
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(1 - e-'^) ± ^(l-e-^)^-8/3Q;'' 



+/'! h-^ \) m 



V{A) can be expressed in a similar way. If A is small, from (12) it follows 
A ~ ift'e*^. Here + (— ) sign corresponds to the case that the r.h.s. in (12) is 
positive (negative), which also corresponds to the case that a' > (a' < 0). 
Then A^ corresponds to cr ^ — oo and V{a) ( V{cr) ) behaves as 

V{a)^V{a)r^-^lnA, • (18) 
a' /x^ 

or 

V{a) = -V{a) ~ -aW . (19) 
On the other hand, when A is large, we find 

Then V{A) is monotonically increasing function for large A if P > and 
approaches to a constant j^^. In order to find the extrema of V{A) and V^(^) 
one differentiate them with respect to A: 

(-^ + 2/5) (A-a' + 2a'\n^) 
V'(A) = V'(A) = . ^ . (21) 
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V'{A) = ii A — a' + 2a' \n ^ = 0, which corresponds to (11). Furthermore 
if a',(3 > there appears a new extremum when —^ + 2/? = or A = 
Aq = \J^- As Aq is not a solution of (11), this may be artificially caused 
by the rescaling. In fact, Aq corresponds to the point that the mapping 
(12) is degenerate, that is, = Q aX A = Aq. Anyway we may discuss 
the (in)stability of the solution in (11) by using the potential V{A). When 
a',/? > 0, let the solution corresponding to (11) be Ai. Combining (18), (20), 
and (21), one finds that if Aq < Ai {Aq > Ai), the solution of (11) (A = Ai) 
is locally stable (instable). On the other hand, if a' < and /3 > 0, there is a 
singularity in V{A) when 1 + ^ + 2/5^4 = 0. It corresponds to cr +oo. Let 
the two solutions (if exist) in (11) be A± {A_ < A^). Combining (18), (20), 
and (21) again, we find that there are three cases: 1. When A2 < A_ < A^, 
A+ (A_) is locally stable (instable). 2. When A^ < A2 < A+, both of A+ 
and A_ are locally stable. 3. When A- < A^ < A2, A- (^+) is locally 
stable (instable). 

When a' > 0, a typical (conceptual) potential is given in Figure 1. There 
may be the following scenario of the inflation and of the present accelerating 
universe: If we start with large curvature A = R = -Rinitiai, the infiation 
occurs due to the large curvature. Since the potential is slowly increasing, 
the curvature rolls slowly down the potential and becomes smaller. If the 
curvature reaches the local minimum A — R — Rq, the curvature will stay 
there with non-trivial value R — Rq ^ 0, which may correspond to the 
present acceleration of the universe. As the potential plays a role of the 
cosmological constant, the ratio w of the effective pressure with respect to 
the effective energy density becomes minus unity, w = —1. As we will show 
later in (46), by fine-tuning the parameters, the cr- field can become heavy 
and decouple. 

3. Cosmic deceleration in Ini? gravity. Let us assume the metric in the 
physical (Jordan) frame (4) is given in the FRW form: 

3 

ds^ = -dt^ + d{tf ^ Qijdx'dx^ . (22) 

In the Einstein frame, the FRW equation looks like 

ml = + \v{a) or 3if| = - \v{a) . (23) 
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Here we distinguish the quantities in the Einstein frame by the subscript E. 
The Hubble parameter He is now defined by He = On the other hand, 
the equation derived by the variation over a is 

= 3 (o- + 3He&) + V'{a) or - 3 (a + 3He&) - V'{(t) . (24) 

We now consider the case that A or the scalar curvature in the Jordan (phys- 
ical) frame is small and the potential is given by (18). Then the solution (in 
the leading order of ig) of the combined equations (23) and (24) is given by 

cr ~ — In — + - In In — •• • , oc t|; + • • • . (25) 
^0 2 to 

The time coordinate in the Euclidean frame is related with the time coor- 

1 1 

dinate t in the (physical) Jordan frame by e^dtE — dt. As a result tQt^ oc t. 
The power law inflation occurs in the physical (Jordan) frame 

2. - 1 

a — e'^ttE OC Ie cx: . (26) 

Since we have a > but a < 0, the deccelerated expansion occurs. In case of 
the original 1/R model in [3], the solution is a oc i^. This might suggest that 
for the model containing i?" (— 1 < n < 0), one may have more moderately 

accelerating universe: d oc t"^ (I < m < 2). 
From (23) and (24), we may evaluate 

a'tl ~ 0{1) . (27) 

Eqs.(26) indicate H = ^ = Since H'^ = ^p, the energy density p corre- 
sponding to a may be deflned as p = ^. Here po is a constant. Denoting 
the pressure of cr by p and substituting the above expressions of H and p to 
the conservation law for the energy-momentum tensor one finds w = ^ — ^, 
which is nothing but that of the radiation. 

One may account for the matter contribution to the energy-momentum 
tensor. When it is dominant compared with the one from a, the obtained 
results are not changed from those in [3] where the possibility of cosmic 
acceleration in 1/R model was established. 

4. Instabilities and corrections to gravitational coupling constant. 

We now discuss the (in)stability of our model under the perturbations. In 
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[7], small gravitational object like the Earth or the Sun in the model [3] is 
considered. It has been shown that the system quickly becomes instable. 

Following to the idea in ref.[15] we start from the action (1), where F is 
given by (10) with m = 2. If the Ricci tensor is not covariantly constant, the 
trace of equation of motion with matter is given by^ 



F(3)(i?)^ , F'(R)R 2F{R) 



K 



2 



°^ + F(2)(i?)^^^^'^ + 3F(2)(i?) 3F(2)(i?) - 6F(2)(i?)^ " ^^^^ 

Here T = T^^. In case of the Einstein gravity, where a' — (3 = 0, the solution 
of Eq.(28) is given by R = Rq = —^T. The perturbation around the solution 
may be addressed^ 

R^R^ + Ri, {\Ri\ < \Ro\) . (29) 
Then by linearizing (28), we obtain 



F(2)(i?o) " " ^ 3F(2)(i?o) 3F(2)(i?o) 3F(2)(/2o) 
+ORi + 2 ^|^)|^°| Vpi?oV^i^i + U{Ro)Ri . (30) 



Here 



F(i)(i?o)F(3)(i?o)i?o F(i)(i?o) _^ 2F(i?o)F(3)(i?o) i?oF(3)(i?o) 



3F(2)(i?o)2 F(2)(i?o) 3F(2)(i?o)2 F(2)(i?o)2 ' 

If (i?o) is negative, the perturbation Ri grows up exponentially with time. 
The system becomes instable. By including the i^^-term, the time for insta- 
bility to occur is significally improved (by the order of 10^*^) [15], compared 

^We should note that the convention of the spacetime signature here is different from 
those in [7]. 

^Duc to the spherical symmetry and the structure of the lagrangian it is enough to 
consider only perturbation of the curvature. Equivalently, one can transform the analysis 
to Jordan (Brans-Dicke) theory. There are no problems with instabilities in the Einstein 
part, while resolution of the instabilities for sigma field which corresponds to curvature 
leads to the same bounds as in below analysis. 
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with the original 1/R model in [3]. If the coefficient of is already fixed by 
some other condition, one can use other HD terms (like i?'^) to eliminate the 
instability completely. 

In (28), we only considered the Ricci scalar. More general equation of 
motion, including all the metric components or Ricci tensor, has the following 
form: 

\g^,F{R) - R^.F{R) - ^^.V^V' {F'{R)) 

+V^V.(F'(i?)) = -^7;. . (32) 

Let the perturbation of the metric be Sgf^u- We may choose a gauge condition 
V^Sg^y = 0. Furthermore we decompose Sg/j^i, into the sum of the trace part 
SG — g^"5g and traceless part dg^^ as 

5g,j.v = + ^guJG . (33) 

For perturbation from the background with constant curvature [R^i, — pgnu) , 
one gets 

SR^R^^ -^SG - W^SG . (34) 

Then if Ri is given, 5G is uniquely determined as in the usual Einstein gravity 
up to the homogeneous part SGh which satisfies —^SGh — V^SGh — 0. For 
the Ricci tensor R^^, we have 

SR,u = - ^ V2%, - ^g^^uSG - ^g^.V'5G - ^V^V,5G . (35) 

Therefore if we use (32), the traceless part Sg^i, can be also uniquely deter- 
mined up to the homogeneous part 6gh^u which satisfies pSgh^u — \^^Sghij.u = 
0, as in the usual Einstein gravity. Since the homogeneous parts SGh and 
^gh^ll' appear in the usual perturbation from the background of the deSitter 
space, they are not related with the (in)stability. Then if SR — Ri is stable, 
whole the metric perturbation 6g^i, is also stable. 

In [8], it has been found that the linearly growing force appears in 
model due to a diffuse source in a locally dcSitter background: 

ds'^^-dt^ + e^'"' J2 {dxf , (36) 

i=l,2,3 
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which is a solution, with a constant curvature Rq = 12H^, of the equa- 
tion corresponding to (28) or (11) for the vacuum case. If we consider 



the perturbation (29) by assuming Ri = Ri{y) {y = e i/y ^j^]^2 3 (^*) )' 
the solution Ri is a sum of two independent solutions fo{y) and f-i{y): 
-Ri(y) — Pifoiv) + P2f-i{y)- First, we consider the case that P ^ a'R^. 
Then U{Ro) = ^ and 

My) ^l-ly + O {y') , U{y) ^-{l-3y + {y')) , (37) 

■J y 

for the vacuum solution. If we assume there is a spherical matter source with 
mass M and the radius Tq as in [8] , the coefficients (3i and (^2 are determined 

by 

p. - ^[fM jr^\ -^{i + ^(yo)}, 

SMCr /:i(l/o)/o(l/o) r' 4MG|/g , 

= -^V-^^y^^ — jrMTS --^^{i+^(^o)}(38) 



Here 16nG = k^. It is assumed that the source exists in y < y^. Since for 
the size of galaxies, y = 10~^ and for the typical distance between galaxies, 
y — 10~^, one may assume yo <ti y <ti I [8]. Then /32 <^ /3i and term with 
P2 may be neglected. If we denote the trace part of the perturbation of the 
metric by h, we find 



(39) 

Then there appears a linear growth as in [8], which might be a phenomeno- 
logical disaster. However, that was the case with large (3. In more general 
case, the equation corresponding to (30) has the following form: 



= 



Ri{y), (40) 



dy"^ y dy Rq 
U{Ro) = — ' ^ ° ■ (41) 
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As clear from Eq.(ll), Rq does not depend on /5. Choosing [3 U(Rq) 
becomes very large and we find i?i — > in the vacuum, which is identical with 
the case of the Einstein gravity with cosmological constant. Then contrary 
to the case of ref.[8], there does not appear the linear growth of h'{y) but 
h'{y) behaves as which docs not conflict with the present cosmology. We 
should also note that the condition 3 7^ is identical with the condition 

that a-field decouples (46), as will be shown below. Hence, modified grav- 
ity with terms growing at small curvature and with higher derivative terms 
important for early time inflation may be viable theory. 

It has been mentioned in ref.[13] that 1/R model which is equivalent to 
some scalar-tensor gravity is ruled out as realistic theory due to the con- 
straints to such theories. As the coupling of a with matter is not small [11], 
we now calculate the square of scalar mass, which is proportional to V"{a). 
We consider the fluctuation from the solution (25). Then 

V"{a) -^-f,. (42) 

E 

1 1 

Here we have used t^t^ oct and (27) are used. The result (42) itself does not 
change from the original 1/R model [3] . Since the Hubble parameter is given 
by if = ^ = ^, in the present universe, we have j ~ ifg- Here Hq ~ 10~^^eV 
is the Hubble parameter of the present universe. Then 

V"ia) tlH', . (43) 

Surely Hq is very small but wc have no restriction (or we have not found it) 
on to- Then if to is very large, the mass of a can be large. Assuming the 
mass is larger than 1 TeV, we have to ~ lO^^eV"^'^ As a't^ ~ 0{1) (27), 
a' ~ (10~^^eV) . This indicates that such class of theories may still pass the 
solar system bounds for scalar-tensor gravity. Moreover, the account of the 

^Thc parameter to should be determined from the initial condition but since there is 
unknown parameter a' , which may differ from (lO~^'^eV) as we have argued, to contains 
the ambiguity coming from a'. As the curvature is small, from (12), we have cr In 
On the other hand from (25), cr^ —In^ — In^ ~ ln(i?oto). As A corresponds to 
the real curvature, one may assume A ~ Hq. Then we obtain a't^ ^ C(l), which is 
identical with (27). In order to determine the value of the parameter a' , one should use 
the information related with the inflation and the matter contents. 
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terms with derivatives of the curvature ^ may permit to pass the solar system 
tests even easier. 

In [13], the PPN (Parametrized Post- Newtonian) parameters have been 
investigated in the Jordan (Brans-Dicke) frame and it has been found that the 
VLBI parameter is above the constraint due the solar system experiments. 
The analysis relied on the mass of the Brans-Dicke scalar, which corresponds 
to a or A in this paper. In the Brans-Dicke form, we find u = 0, that is, 
there is no kinetic term for a. In the current limit, u > 3500 [14]. Then the 
Brans-Dicke scalar should be heavy to avoid the problem. In the Minkowski 
background, it has been found that the mass is too light and unnatural. In 
our case, the action correponding to (10) can be rewritten as 



5 ^ J_ / d'x^ 



, (44) 



which is also the Brans-Dicke form with a; = 0. We should note, however, 
the present universe is not exactly Minkowski but accelerating. The mass of 
the BD scalar is given by (43). Then as discussed, the mass is determined by 
the parameter of the integration to? which may be determined dynamically. 
Such a mass may be quite large and then such theory does not violate solar 
system test as in [13]. 

One may consider the case that the present universe corresponds to the 
solution (11). In such a case, by tuning the parameters /? and a', the mass 
of a can be made large again. Let us write the solution(s) of (11) as Ai for 
a' > case and A± for a' < case. Then 



da^ 



A=Ai,A± 



'da\ d^V{A) \ 
dAj dA^ J 



A=Ai,A± 



{l + ^ + 2/3A) {-^ + 2/3) 



A=Ai,A± 



Choosing 



a 



2^2 



(45) 



(46) 



A=Ai,A± 



^The correspondent scalar-tensor theory becomes higher derivative one as it follows 
from second section. 
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for a' > or 

/^-^ - 77T + 



2A 2^2 



(47) 

A=AuA± 



for a' < , the mass of a becomes large again. Thus, HD term may help to 
pass the solar system tests for In i? (or 1/ R) gravity. 

5. Discussion. Finally, we calculate the corrections to gravitational cou- 
pling constant. The easiest way is to consider the perturbation around the 
constant curvature solution (11). We now write the metric as g'^l^} corre- 
sponding the solution and the constant scalar curvature as R^^'' — ||. The 

metric is splitted into the background part gf} and the perturbation h^i,: 

9^^^ = gtu + V ■ (48) 
The following gauge conditions are imposed: 

= ^(0)'^'' V = V^'^)^ V ■ (49) 

The first condition is chosen to simplify our discussion as graviton is spin 2 
field. Then 



(50) 



Using the equation = R^'^^ + 2a' In ^ — a', the expanded action 5* has the 
following form: 



f^) / d^.v^ (^ - ijA^^ft- + iA-D'")/., 



12 P I I " " \P 2^ "" 4 

(51) 

Here the total derivative term is dropped. In case of the Einstein action with 
a positive cosmological constant 

^^ = ^/v^(^-^) > (52) 
which has deSitter solution, the corresponding linearized action is: 

~SE = lr,j d'x^/^ (^^ - ^^K^h^^ + ^/I'^'^D^^) v) ■ (53) 



13 



Comparing (53) with (51), the gravitational constant k is renormahzed as 



1 . 1 f,.»''\^P) . (54) 



12 P 

The Newton potential will be modified respectively. When Rq = ^ cor- 
responds to the rate of the expansion of the present universe, one gets 
-Ro ~ ^ ~ A*^ ~ lO^^^eV. In (54), /3-dependent term is of the same or- 
der term as second term when scalar has large mass. With the assumption 
a' ~ (10~^'^eV) the correction may be significant. We may, however, con- 
sider the case a' is much smaller than 10~^^eV, since there seems to be no 
constraint for the value of a' itself. We only have a constraint (27) and to 
can be very large. In such the correction to the Newton constant can 

be very small and the starting theory has acceptable newtonian limit. 
One may discuss further generalizations of modified gravity like 

F(yl) = A + hn— J . (55) 

Here we restrict n by n > — 1 (m is arbitrary) in order that the second term 
^ could be more dominant than the Einstein term when A is small, n and m 

can be fractional or irrational numbers in general. 

When the physical scalar curvature A = R is small, we find 

/ A\™ 1 / a \^ 

e-<^ 7n^-"-^ [In—] , A^ {--fne")^ 



y(a) - ( 1 + - ) (-7n)^e^'' f - ) , (56) 



when n ^ and 



A* / \ n -|- 1 

1 n+2 / a 



V{a) ~ 76^^" , (57) 



^The sum of such terms (whose coefficients could be constrained by the condition of 
avoiding the Hnear growing gravitational force) may be considered. The presence of E? at 
large curvature is supposed. 
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when n = 0. With the similar procedure as in the previous sections, 



Ti ~ C(l) . (58) 



7 
"•0 



As a result 



3(»+l) 



2 

32- (n+l)(2n+l) 



t ~ , aE ~ , a ~ t — . (59) 



This does not depend on m. The logarithmic factor is almost irrelevant. We 
also find the efi^ective w for the cr-field is 

6//' + Iv — 1 , . 

" = - 3(n + l)(2n + l) - ^'^^ 

Then w can be negative if 

1 —7 + \/73 
-l<n<--orn> — = 0.1287---. (61) 

Prom (59), the condition that the universe could accelerate is > 1, 

that is: 

n > ^ ^ = 0.366 • • • . (62) 

Clearly, the effective dark energy w may be within the existing bounds. 

Thus, we demonstrated that modified gravity with Ini? or i?~"(lni?)"* 
terms may be responsible for the current acceleration of the universe. Hence, 
like the simplest 1/ R modified gravity this provides the gravitational alterna- 
tive for dark energy. Moreover, the presence of HD terms like (which may 
be responsible for early time inflation) helps to pass the existing arguments 
(instabilities, solar system tests) against such modification of the Einstein 
gravity. The theory may also have the wcU-acccptable newtonian limit. It 
is clear that much more work is requiered to (dis)prove that one of the ver- 
sions of such modified gravity is currently realistic theory. Nevertheless, the 
fine-tuning of parameters of modified gravity to provide the effective gravi- 
tational dark energy looks more promising than the introduction by hands 
some mysterious fiuid. 
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'initial 



Figure 1: A typical potential when a' > 0. We may start with large curvature 
A = R = i?initiai (inflation). Then the curvature rolls down the potential 
slowly and stops at the small curvature A — R — Rq (the present accelerating 
unvierse) . 
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